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MULTI-CAST COMMUNICATION SYSTEM 
AND METHOD OF ESTIMATING CHANNEL IMPULSE RESPONSES 

THEREIN 

5 Background to the Invention 

This invention relates, in general, to a multi-cast communication system and a 
method of estimating channel impulse response (IR) therein, and is particularly, 
but not exclusively, applicable to communication environments employing space- 
time coding. The present invention is also applicable, without imposing limitation, 
10 to code division multiple access (CDMA) schemes, orthogonal frequency division 
multiplexing (OFDM) or the global system for mobile (GSM) communication, and 
also to systems having base stations configured to transmit on either a time- 
aligned or unsynchronised basis. 

15 Summary of the Prior Art 

GSM, CDMA, and OFDM systems and those systems using space-time coding 
usually require estimation of impulse responses (IR's) between a terminal 
antenna and several base antennas, especially in instances when the relative 
signal strengths of the base stations are similar. As will be understood, the 

20 physical channel through which propagation occurs can have a severely 
detrimental effect on the ability to recover data accurately, especially with 
increasing data rates. Consequently, by knowing the IR, compensation for 
channel-induced defects can be applied to improve accuracy in data recovery, 
e.g. through interference cancellation techniques. More especially, with multiple 

25 transmit elements at a base station, irrespective of whether there is more than 
one serving base station, multipath interference results in an inability to resolve 
individual channels. Indeed, even orthogonally structured data (such as different 
and time dispersed training sequences) can become cross-correlated (in the 
face of multipath interference) and hence unusable to resolve individual 

30 channels/paths in the context of a deterministic IR approach. Furthermore, it will 
be appreciated that IR determination is particularly taxing at a handoff point at 
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cell boundaries of, potentially, two or three serving base stations (or Node Bs in 
the context of CDMA and the Universal Mobile Telecommunications System 
(UMTS)). At cell-boundary handover, the issue of I R is complicated by that fact 
that downlink carrier to interference ratios (CIRs) can be zero or negative. 

When available, impulse responses can also be used to assist in downlink soft- 
handover. 

The use of matched filtering of standard training sequences to determine 
0 channel IR's is considered adequate to drive fingers of a RAKE filter in CDMA 
applications, such as IS-95 and CDMA2000, given that incoherent soft handoff is 
also allowed. However, more recent cellular system proposals will stretch 
matched filtering techniques to their limits and beyond due to the use of higher 
data throughputs (with more bits per symbol, e.g. higher quadrature amplitude 
5 modulation QAM levels) requiring higher energy per bit/noise floor (E b /N 0 ) than 
existing systems, such as IS-95. Such higher data rate cellular systems include 
the Enhanced Data-rate for GSM Evolution (EDGE) system. 

In any event, interference cancellation, coherent soft hand-off or space-time 
20 coding generally requires a more accurate estimate of channel propagation 
conditions (than that required with RAKE filtering or matched filtering). In this 
respect, a long history of interference cancelling base station algorithms (that 
yield indifferent null-steering performance) lends support to this view. 

25 In terms of proposed 3G (third generation) cellular systems, it is known to use 
training sequences with a cyclic prefix for pilot tones in orthogonal frequency 
division multiplexing (OFDM) and CDMA. When a channel has an impulse 
response of known maximum duration, the use of Steiner codes allows re-use of 
the same training sequence with a time-offset method to measure 

30 simultaneously the channel impulse responses of several users or several 
antenna with a single fast Fourier transform (FFT) operation Steiner codes are 
discussed in the paper "A comparison of uplink channel estimation techniques 
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for MC/JD-CDMA transmission systems" by B. Steiner and R. Valentin, 
Proceedings IEEE 5th International Symposium on Spread Spectrum 
Techniques and Applications 1998, Volume: 2, pp. 640-646. Indeed, cyclic 
Steiner codes are found in the pilot tones for the time division duplex (TDD) 

5 component of European UMTS, as presented in the ETSI 3GPP Document 
TS25.211 "Transport channels and physical channels" v.3.2.0, 1999. In terms of 
a frame structure for a pilot tone, random data sectors sandwich an observation 
window that is preceded by a cyclic prefix of a latter portion of the observation 
window. Provided that the duration of the cyclic prefix exceeds the channel 

10 impulse response, then data in the observation window is a function purely of 
the pilot sequence (comprised from the combination of the cyclic prefix and data 
in the observation window) and is not corrupted by spurious unknown data in 
any way. 

15 The 3G Partnership Project (3GPP) have defined training sequences for TDD 
(time division duplex) having overall chip lengths of 256 and 512 chips; the 
overall lengths are derived from 192 chips and 464 chips of basic pseudo- 
random number sequence codes with, respectively, an additional cyclic prefix of 
either 64 chips and 114 chips. Each training sequences is designed to allow 

20 channel impulse response estimation for differing numbers of users, namely 
three users in the case of the 256-chip training sequence and eight users in the 
case of the 512-chip training sequence. 

The training sequences are utilised in the receiving unit to estimate the channel 
25 impulse response based on a complex cross-correlation (in real and imaginary 
phase and amplitude components) between received chips and a local replica of 
the training sequence. In this regard, it is usual, on a per channel basis, to take a 
correlation of the channel impulse response (h) with the sequence (s) from the 
transmitter. More particularly, from a single base station having multiple transmit 
30 elements, cyclic offsetting of Steiner codes allows utilisation of a fast Fourier 
transform (FFT) technique to solve individual channel impulse responses. 
Steiner cyclic pilot codes can therefore be used in estimating, with a single 
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correlator, channel impulse responses of multiple users that do not mutually 
interfere. Steiner codes may be Gold codes. 

As will be understood, the requirement for cyclic redundancy to mitigate data 
5 corruption caused by multipath may be obviated provided that the length of the 
transmitted training sequence (or pilot tone) is sufficiently long. In generality, 
cyclic redundancy may be avoided if a training sequence contains a sufficient 
number of chips, with the overall length of the training sequence determined by 
the prevailing dispersion conditions associated with the channel. Additionally, for 
10 estimating CIR, it is also necessary to consider the physical state of the 
receiving unit, since the channel for a slow moving or stationary receiving unit 
(e.g. a handheld device) is ostensibly stable, i.e. constant. 

By way of practical explanation of Steiner code operation from a dual element 
15 array of a base station transmitter, it will be appreciated that, in estimating the 
two channels from the two radiating elements to a single receive antenna, the 
two radiating elements use a common generic pseudorandom number (PN) 
training sequence {Si, §2} which is end-around shifted (in this case by half a 
block) to {§2, si} for the second element. The two channel impulse responses are 
20 {£2} and the receiver antenna sees the superposition of the two convolved 

sequences: 

{y} = ,s 2 ) ® i\ ) + {si> s i )®{ } h} 
where <8> means a discrete-time convolution operation. When the receiver output 
is correlated with the PN sequence of the first antenna, the impulse responses of 

25 the two distinct channels are recovered in superposed form with a relative time 
shift of half a block. If the maximum probable duration of the impulse responses 
is finite and known and less than half the PN code length, the two channel 
estimates are non-overlapping and orthogonal and can be recovered by slicing 
into two halves an output of a correlator block. Clearly the process extends 

30 easily to estimating M channels by splitting the generic PN sequence into M 
segments {§1, §2,.. sm}, where {s} means a short symbol sequence, and wherein 
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transmissions from the base antenna transmit elements follow the sequences 
{§1, Sm}, {S/W, Si,.,. Sw-i}i--.{Sw-i* 

At the receive antenna, the processing block is correlated with only the first 
5 copy, {si, s 2 ,.. sm} of the PN sequence whereupon, at the correlator output, the 
M different channel impulse responses separate out in time and appear in 
sequence. Clearly, the cyclic head is different for each downlink path, but this is 
purely a transmitter burden, whereas the receiver samples the waveform and 
continues to see the same rotated training sequence. In this way, a common 
10 discrete Fourier Transform operation can be used to resolve multiple downlink 
channel impulse responses. 

As previously indicated, in order to avoid contamination with unknown data 
which usually immediately precedes the training sequence, a cyclic copy of the 
15 end of each rotated sequence may be prefixed at the start of the pilot block. In 
other words, when working with, say, 512-point processing blocks and when the 
duration of the channel impulse response is K samples, the total extent of the 
transmitted pilot burst is 512+K samples of which the first K contaminated 
samples are discarded at the receiver. 

20 

The maximum number M of channel impulse responses that can be estimated 
by the Steiner method is; 

Ml" 

K 

where N is the length of the training sequence (samples) and K is the maximum 
25 likely number of samples in the channel impulse response. 

Unfortunately, with Steiner codes, the resolution of base stations in different 
cells is much less attractive. Specifically, whilst it would be possible to continue 
to expand the Steiner method with ever increasing numbers of built in channels 
30 M (e.g. perhaps twice or three times the number needed for one base station), 
estimating more channels needs longer training sequences and synchronisation 
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of the relevant base stations. Increasing training sequence length, however, is 
not the main difficulty; a bigger problem is that this Steiner solution would 
effectively constrain all the bases to use one and the same fundamental PN 
sequence and this would not allow addressed units to discriminate between 
bases by different codes when they are searching for hand-off possibilities. 
Steiner is also inflexible since there is contention over the order in which the 
bases should cycle their training sequences. This contention will appear at all six 
interfaces of the base's cell with other cells in a hexagonal layout and so 
sequence cycling must be set up to be compatible with all adjacent sites that 
could be affected. 

As will now be appreciated, communication systems demand the resolution of 
multiple channels that combine to provide a composite channel impulse 
response C from observed training sequences S transmitted by several base 
stations each with multiple antenna elements. In this regard, a time domain 
channel estimation hi(t) for pilot sequence stf) convolved through channel h(t) 
that is subject to noise n(t) may be obtained through a time-reversed correlation 
of the received sequence with a matched filter having the form Si*(-t). In the 
frequency domain, the channel estimation takes the form Hi(g>^. 

For the purposes of channel estimation, matched filtering of the PN sequence in 
the receiver can be performed through fast Fourier Transform (FFT) techniques. 
If the received data block of N samples is {yi, y 2 ... yn} (which could be a 
sequence like {si, s 2 ... s M }), then FFT operation yields: 

FFT 
FFT 

In other words, matched filtering operation in the time domain is equivalent to a 
point-by-point multiplication in the discrete Fourier domain, namely: 

H^=Y k -Sl k = 0...2"-\ 
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with two different channel estimates contained in the semi-sequences of {/7>, 
namely: 

{h)i = A/2— Vi) 

Of course, once in the frequency domain, sidelobes present in the frequency 
5 spectrum in a channel estimation circuit containing a matched filter may be 
minimised through the use of a Wiener (least squares) filter technique. 



The Wiener least squares filter provides an equalisation technique that utilises a 
modified inverse filter that controls the white noise response of the filter, i.e. the 

10 undesired enhancement of thermal noise from the antenna. The paper "Smart 
Antennas for Third Generation Mobile Radio Systems" by Martin Haardt 
(Siemens), Stanford Colloquium on Smart Antennas", July 1999, describes 
channel equalisation in terms of the Wiener filtering response, and explores how 
two or more receiving antennas can be included in the description of the 

1 5 received data. A further paper by H Sari et al titled "Transmission Techniques for 
Digital TV Broadcasting", IEEE communications magazine 33(2) Feb. 1995, 
discusses channel equalisation in the context of the Wiener filter mechanism. 
Space-time transmit diversity (STTD) is included in the standards for 3G cellular 
systems in the European UMTS Terrestrial Radio Access (UTRA) system and 

20 similar systems in the American CDMA2000 proposals. Consequently, it is 
important for these STTD systems to implement efficiently data recovery and 
equalisation algorithms. Utilisation of the Wiener filter in a communication 
system environment is further discussed in the Applicant's co-pending European 
Patent application 01300520.2 (Applicant's reference 11872ID-Hudson) having 

25 both a priority claim to US patent application serial number 09/488,721 and 
corresponding to the continuation-in-part application (Applicant's reference 
11872IDUS04I-Hudson) filed in respect thereof, all incorporated herein by 
reference. 
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Equalisers and channel whiteners are generally not well conditioned, especially 
if the channel has zeros or deep minima in its frequency response. The stable 
minimum mean square error (MMSE) Wiener filter solution for the channel is 
therefore: 

y ■ 9* 



where a 2 is the variance of the thermal noise level in the frequency domain. If 
the DFT is orthogonal, i.e. a unitary matrix operation, then this a 2 value is the 
same as the time domain noise variance per sample, but most FFT algorithms 
apply some form of scaling for which allowance must be made. 

10 

Use of the Wiener filter solution results in residual errors caused mainly mainly 
by thermal noise in the sidelobes (and not error due to imperfections in the code 
sequence); this can be contrasted with PN sequence autocorrelation sidelobes 
experienced in a matched filter environment. 

15 

The paper by J. Blantz et al on the "Performance of a cellular hybrid C/TDMA 
mobile radio system applying joint detection and coherent receiver antenna 
diversity", IEEE J. Selected Areas in Comms. 12(4), May 1994, pp.568-579, 
describes a multi-user detection algorithm. A second paper by A. Klein et al on 
20 "Zero forcing and minimum mean square error equalisation for multiuser 
detection in code-division multiple access channels", IEEE Trans. Veh. Tech. 
45(2), May 1996, pp.276-287, provides further context to the present invention. 

In overview, therefore, PN training sequences sent from multiple base stations 
25 (or Node Bs) employing one or more transmit elements are subject to multipath 
that results in code cross-correlation and an inability at a receiver to resolve the 
individual channels and establish the individual channel-specific impulse 
responses. More specifically, whilst time alignment of base station transmissions 
may produce a summation of signals at the receiver (when employing 
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appropriate windowing on a chip-by-chip basis with respect to identifiable 
correlation spikes), the receiver is only able to detect a composite channel 
impulse response that is unlikely to reflect accurately any of the actual 
transmission paths. Moreover, whilst the composite channel impulse response 

5 may be sufficient in the context of soft handover (in IS-95, for example), the 
composite channel impulse response is generally insufficient in third generation 
systems, including systems offering space-time coding. Furthermore, there is a 
reticence shown by service providers to provide synchronicity between base 
station transmissions since synchronised transmissions increase infrastructure 

10 costs, such as through the necessary provision of an accurate timing reference. 

Summary of the Invention 

According to a first aspect of the invention there is provided a method of 
determining channel impulse responses of a plurality of channels to a 

15 communication device, the method comprising: performing transform operations 
on both a replica of a signal sequence s n and a received training sequence y n by 
the communication device in at least one burst, the received training sequence 
y n being the signal sequence as received through a channel, the transform 
operations arranged to generate a multiplicity of signal sequence frequency bins 

20 and a multiplicity of training sequence frequency bins; performing point-by-point 
operations between corresponding signal sequence frequency bins and training 
sequence frequency bins; and concatenating the point-by-point operations 
associated with the channel to provide a composite frequency response for the 
channel, the composite frequency response allowing, in the time domain, 

25 generation of the channel impulse response for the channel. 

It is preferred to use a set of matrix operations in the frequency domain to 
resolve channels to the communication device from multiple transmitting 
stations, the matrix operations providing solvable linear equations for the 
30 channel impulse response H and the signal sequence S expressible in a matrix- 
vector form at each frequency bin. 
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A minimum mean square error (MMSE) estimation technique is preferable to 
solve the linear equations, although other mathematical techniques are possible 
(as will be appreciated). 

The MMSE estimation technique preferably employs a Weiner filtering operation 
providing: 

H = (s H S + Cov(B)Ys H Y 

where: Cov(B) is a covariance matrix of thermal noise; S H Y is a matched filter 
operation arranged to minimize the signal to noise ratio but not to remove mutual 
interference between channel estimates; and (S H S+Cov(B))'' i is a decoupling 
matric that removes the coupling between different channels caused by any non- 
ideal nature of the training sequences. 

In a preferred embodiment, the number of bursts sent to the communication 
device from each transmitting unit in communication contact therewith is 
calculated as a multiplication of: a number of transmitting elements in a transmit 
array of a transmitting unit; and a number of transmitting units in communication 
contact with the communication device. 

In a second aspect of the present invention there is provided a method of 
determining channel impulse responses of channels incident to a communication 
device, the method comprising: and transmitting multiple quasi-orthogonal 
pseudo-noise sequences as bursts from multiple base stations each having at 
least one transmit element, successive bursts providing an extended training 
sequence for use in channel estimation at the communication device; applying a 
Wiener frequency domain MMSE deconvolution with frequency domain spatial 
decoupling matrices to generate channel impulse response estimates for the 
channels. 

In another aspect of the present invention there is provided a method of 
determining channel impulse responses of a plurality of channels established 
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between a plurality of transmitting elements and a communication device in a 
communication system, the method comprising: substantially simultaneously 
transmitting different training bursts from, each of the plurality of transmitting 
elements, each burst having a length at least as long as a maximum channel 

5 duration in the communication system multiplied by a number corresponding to 
the plurality of transmitting elements; recovering at the communication device a 
signal sequence y n from the different training bursts s n ] and resolving the 
plurality of channels to recover associated channel impulse responses H for 
each channel by solving an algebraic matrix operation expressed in matrix- 

10 vector form as Y=SH r where: S is a matrix of partial training bursts for each 
channel, each training burst segmented into N pieces in the time domain; Y is a 
vector of a received signal sequence; and H is a concatenation of different 
channel impulse response vectors. 

15 A further aspect of the present invention provides a computer program product 
for a processor within a receiver dfevice, the computer program product 
comprising: code that performs transform operations on both a replica of a signal 
sequence s n and a received training sequence y n received by the communication 
device in at least one burst, the received training sequence y n being the signal 

20 sequence as received through a channel, the transform operations arranged to 
generate a multiplicity of signal sequence frequency bins and a multiplicity of 
training sequence frequency bins; code that performs point-by-point operations 
between corresponding signal sequence frequency bins and training sequence 
frequency bins; and code that concatenates the point-by-point operations 

25 associated with the channel to provide a composite frequency response for the 
channel, the composite frequency response allowing, in the time domain, 
generation of the channel impulse response for the channel; wherein the codes 
reside in a computer readable medium. 

30 According to another aspect of the present invention there is provided a 
communication device having a receiver coupled, in use, to receive a plurality of 
channels supporting a signal sequence y n and training sequence bursts, the 
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communication device having: a signal processing platform to perform transform 
operations on both a replica of a signal sequence s n and a received training 
sequence y„ received by the communication device in at least one burst, the 
received training sequence y„ being the signal sequence as received through a 

5 channel, the transform operations arranged to generate a multiplicity of signal 
sequence frequency bins and a multiplicity of training sequence frequency bins; 
the signal processing platform arranged to perform point-by-point operations 
between corresponding signal sequence frequency bins and training sequence 
frequency bins; and the signal processing platform further arranged to 

10 concatenate the point-by-point operations associated with the channel to provide 
a composite frequency response for the channel, the composite frequency 
response allowing, in the time domain, generation of the channel impulse 
response for the channel. 

15 Yet another aspect of the present invention provides a communication receiver 
comprising: means for receiving, in use, multiple quasi-orthogonal pseudo-noise 
sequences as bursts from multiple base stations each having at least one 
transmit element, successive bursts providing an extended training sequence for 
use in channel estimation at the communication receiver; and means for 

20 applying a Wiener frequency domain MMSE deconvolution with frequency 
domain spatial decoupling matrices to generate channel impulse response 
estimates for the channels. 

In one aspect of the present invention there is provided a communication device 
25 operational to receive a plurality of training sequences on a plurality of channels 
and a signal sequence y„, the communication device comprising: a receiver for 
substantially simultaneously receiving, in use, different training bursts from each 
of the plurality of channels emanating from a plurality of transmit elements, each 
burst having a length at least as long as a maximum channel duration multiplied 
30 by a number corresponding to the plurality of transmit elements; recovery 
circuitry for recovering, in use, the signal sequence y n and the different training 
bursts s„; and a processor arranged to resolve the plurality of channels to 
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recover associated channel impulse responses Hfor each channel by solving an 
algebraic matrix operation expressed in matrix-vector form as Y=SH, where: S is 
a matrix of partial training bursts for each channel, each training burst 
segmented into N pieces in the time domain; Y is a vector of a received signal 
5 sequence; and H is a concatenation of different channel impulse response 
vectors. 

According to the present invention, a base station of a communication system 
may be provided in which base station a transmitter chain is arranged to transmit 

10 multiple quasi-orthogonal pseudo-noise sequences as training bursts s„ from at 
least one transmit element and further arranged to transmit a signal sequence 
y n , successive training bursts providing an extended training sequence for use in 
channel estimation at a communication device of the communication system, the 
transmitter chain substantially simultaneously transmitting, in use, different 

1 5 training bursts from each of the at least one transmit element, each training burst 
having a length at least as long as a maximum channel duration in the 
communication system multiplied by a number corresponding to a plurality of 
channels to the communication device, the extended training sequence and the 
signal sequence y n providing a resolution mechanism to the communication 

20 device allowing the communication device to resolving the plurality of channels 
to recover associated channel impulse responses H for each channel by solving 
an algebraic matrix operation expressed in matrix-vector form as Y=SH, where: 
S is a matrix of partial training bursts for each channel, each training burst 
segmented into N pieces in the time domain; Y is a vector of a received signal 

25 sequence; and H is a concatenation of different channel impulse response 
vectors. 

Advantageously, the present invention provides a mechanism for resolving a 
composite channel into a plurality of sub-channels with which specific individual 
30 channel impulse responses can be associated with and calculated using a 
numerically efficient fast transform technique. In other words, irrespective of 
whether or not there are synchronized base station transmission to an 
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addressed unit from multiple base stations having one or more radiating 
element individual channels can be resolved to generate channel impulse 
responses to aid in data recovery. Beneficially, the present invention is 
particularly applicable in the estimation of the channel impulse responses 
between multi-element base stations (in different cells) and the receive antenna 
array of a subscriber unit located at a 2-way or 3-way handoff point where 
interference is generally very high. 

It can be shown that for a dispersive channel having two transmitting base 
stations, the present invention provides increased accuracy in channel 
estimation, with smaller errors experienced than in an uncoupled case. Indeed, 
tests have shown that normalised mobile station errors in channel estimation for 
the present invention can achieve results of 14.8 decibels (compared with 6.1dB 
and 7.1dB for techniques solely employing a matched filter and Weiner filter, 
respectively, in channel impulse response estimation). The present invention, 
even in conditions of varying base station transmit power, provides an improved 
accuracy in channel impulse response estimation, which accuracy is reflected in 
a general measured improvement of about 7dB over a range of transmit signal 
strengths of ±20dB for second base station (relative to a carrier signal from a 
first base station being set at +3dB above thermal noise at the clip level). The 
present invention is therefore able to operate in interference and soft handover 
environments where interference becomes a problem. 

A solution for estimating several downlink channels from different base antennas 
to a single terminal antenna by the use of multiple training bursts has been 
described. The method involves the use of NxN decoupling matrices in the 
Fourier domain to solve for N channels and gives superior results compared with 
attempts to use matched filters designed for added white Gaussian noise 
(AWGN) channels. Interference cancellation, space-time coding, and coherent 
soft hand-off can be enhanced by the use of the present invention. 

Brief Description of the Drawings 
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At least one exemplary embodiment of the present invention will now be 
described with reference to the following drawings, in which: 

FIG. 1 is a schematic block diagram of a basic transmitter and receiver 

architecture 10 in which processing circuitry may be adapted to support the 
5 underlying principles of the present invention; 

FIG. 2 is a diagrammatic illustration of a typical cellular communication 

system; 

FIG. 3 shows a mechanism and schematic block diagram for channel 
impulse response estimation that may be adapted to support the principles 
10 of the preferred embodiments of the present invention; 

FIG. 4 is a transmission regime for a training sequence according to a 
preferred embodiment of the present invention; 

FIG. 5 a schematic block diagram of the channel estimation process and 
architecture according to a preferred embodiment of the present invention; 
15 and 

FIG. 6 is a Steiner code variant of FIG. 4. 



Detailed Description of a Preferred Embodiment 

Before discussing the preferred embodiment of the present invention, it should 
20 be noted that consideration has been given to employing Steiner codes in 
combination with Wiener MMSE estimation in the Fourier domain to obtain a low 
cost multi-user channel estimation method that is highly efficient for 
simultaneously estimating a number of channels on the downlink for CDMA 
systems. Steiner codes may also be used in OFDM systems though, in this 
25 case, it may be preferred to use the intrinsic OFDM symbol for training 
sequences. Given only prior knowledge of the maximum extent of the channel 
impulse response, the same basic pilot sequence can be time shifted and re- 
used (with properly determined cyclic headers) to function for a number of 
different downlink channels, and only one FFT correlation process is needed to 
30 solve for all the channels simultaneously. The technique is ideal for wideband fat 
pipe systems where the same user signal is radiated from different antennas of a 
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transmit diversity antenna or a space-time coded antenna system. Channel 
estimation by this method is accurate enough to allow adequate interference 
cancellation of intercell interference, which is usually quite difficult to achieve in 
CDMA systems due to low grade estimates of the interference parameters. 

Referring now to FIG. 1, there is shown a basic transmitter and receiver 
architecture 10 in which processing logic or circuitry may be adapted to support 
the underlying principles of the present invention. Transmit and receive chains 
100, 101 may be combined in a transceiver arrangement of a station, such as a 
base station/Node B or a mobile device of a cellular radio system or the like. 
Individual transmit and receive chains could, of course, be implemented in 
isolation, whereby unidirectional (as opposed to bi-directional) communication is 
supported across, for example, an air-interface, such as a CDMA or 3G radio 
environment (although the present invention is more widely applicable to multi- 
point transmission systems). 

For the sake of explanation only, the following description of the preferred 
embodiment will be based on a wideband CDMA system and so references to 
"CDMA" should be construed broadly and without specific limitation. 

In transmit chain 100, signal source 102 generates signals 104 which are 
encoded by encoder 106, such as a half rate turbo coder. Encoded signals 107 
are modulated and subjected to CDMA spreading by modulator/spreader 108 
and then filtered by pulse shaper digital filter 110. Following filtering, a digital to 
analog (D/A) converter is arranged to produce analog signals 114 which are 
generally filtered in a second filter 1 16 to limit any aliasing effect at transmission. 
A mixer 118 up-converts baseband signals 120 to transmission frequencies prior 
to amplification in power amplifier 122 and transmission from antenna 124. Of 
course, as will be appreciated, some of the functional blocks in the transmit 
chain can be realised as discrete or collocated into a single chip or processing 
block. 
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Transmitted signals 126 are received by receive antenna 128 of an addressed 
station. Received signals 130 are communicated to a down-converting mixer 132 
via an amplifier 134, which down-converting mixer typically produces an 
intermediate frequency signal for initial signal processing. A filter 134 isolates a 

5 desired frequency spectrum 136 for signal processing, which frequency 
spectrum is converted to a digital domain representation 137 by analog to digital 
(A-D) converter 138. The digital domain representation 137 is digitally filtered in 
digital channel pulse shaper filter 140 and then applied in series to a combined 
demodulator/equaliser slicer and RAKE despreader 142 and then a channel 

10 decoder 144. Recovered information 145 is ultimately received by terminal 
receiving equipment (TRE) 146. It will be understood that channel filtering 
generally distorts the channel impulse response and the channel filter becomes 
part of the unknown channel. 

15 The transmit and receive chains 100, 101 are invariably micro-processor (uP) 
controlled operations. By way of illustration, a controller 150 and associated 
memory 152 (which may be a combination of random access (RAM) and read 
only (ROM)) is shown as being common to both the transmit and receive chains 
100, 101. The memory 152 stores control algorithms, such as equalisation and 

20 modulation algorithms, and also operands and incident data is subject to signal 
processing requirements. Moreover, common circuitry between the transmit and 
receive chains can be shared to reduce component count, subject to the 
inclusion of suitable switching and regulating circuitry. Additionally, as will be 
understood, at least some of the functions within one or both of the transmit and 

25 receive chains 100, 101 can be implemented either entirely or partially as code, 
e.g. equalisation and demodulation. It is these code blocks, application specific 
integrated circuits or memory modules that can be adapted or up-graded to 
support the present invention. 

30 Wideband CDMA-specific components can be replaced or substituted to support 
an arbitrarily selected communication protocol, e.g. an FDD environment 
employing TDM techniques or a multicarrier environment. 
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FIG. 2 is a diagrammatic illustration of a typical cellular communication system 
170 showing a subscriber unit 172 (preferably having multiple receive elements 
174-180) at cell boundaries between three potentially serving base stations 182- 

5 186. For an equidistant displacement of the subscriber unit from the base 
stations, one could expect that transmit signal strengths to be roughly equal, 
subject to the respective paths and the associated physical environments. The 
general purpose of FIG. 2 is to illustrate a physical environment in which the 
present invention finds exemplary application. For purposes of clarity, the 

10 cellular communication system 170 is shown with three cells 188-192 and a 
single subscriber unit 172, but clearly cellular systems may include many cells in 
multiple layers, with each cell servicing a multiplicity of subscriber units. At the 
hand-off point midway between two or three base stations, it is necessary to 
distinguish the channel impulse responses for a doubled or tripled number of 

15 channels from multiple base antenna elements. If the CIRs are not estimated 
jointly, the carrier to interference (C:l) ratio is very low within the training 
sequences and the channel estimates will be correspondingly poor. Given that 
each base station always has several TX element channels to estimate it is 
reasonable to use the Steiner method in the intracell mode, allocating one PN 

20 sequence per base station with cyclic rotation to distinguish between the bases' 
antenna elements, as described previously. 

Incoming signals 126 to the subscriber unit 172 are received from the antenna 
128 at a transmission rate of 1 /Tch of the system. By way of practical 

25 implementation, an entire slot of a frame is sampled (A/D converted) at the 
Nyquist rate and read into memory 152 (of FIG. 1) prior to commencement of 
information recovery. The incoming signals 126 can be received from any digital 
source of signals and are not restricted to wireless communications, with the 
base stations utilising single or multiple radiating elements in each transmit 

30 antenna array. Signal processing of incoming signals 126 for optimised data 
recovery is clearly subject to knowledge of the channel impulse response (i.e. 
the "channel spectrum"), with this determined from training sequences sent from 
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the base stations according to a predetermined regime, e.g. once per slot, once 
per frame or on any other appropriate basis. An exemplary mechanism for data 
recovery is described in the aforementioned European and US patent 
applications (e.g. EP 01300520.2). 

As regards the channel spectrum, reference is made to FIG. 3 that shows a 
mechanism for channel impulse response estimation that is adapted to support 
the principles of the preferred embodiment of the present invention. A training 
sequence 250, preferably located in and extractable from each slot of a frame, 
falls incident on receive antenna 128, which training sequence is demodulated in 
the receive chain (as will be readily appreciated) to recover the random 
sequence y k , e.g. a generated pseudo-random number sequence. The form and 
structure of the training sequence will be described in detail later, since the 
structure is inexorably associated with an ability to resolve individual channels 
for channel impulse response estimation according to the present invention. 

In overview and without initially discussing the preferred embodiment, the 
recovered random sequence y k is subjected to a DFT function 252 to produce a 
channel sequence spectrum Y k 254. Within the memory 152 associated with the 
receive chain there is stored a replica of the original training sequence s k 260, 
which original training sequence s k 260 has a corresponding known sequence 
spectrum S k 262 obtained by subjecting the original training sequence s k 260 to 
an FFT function. The known sequence spectrum S k 262 may simply be stored in 
the memory since it is this sequence spectrum S k 262 that is used to assess the 
channel impulse response. 

The DFT function 252 may be performed by a specific digital signal processor 
(DSP) or the processor 150 of FIG. 1. The channel spectrum therefore assumes 
the form of an x-point, e.g. 512-point, DFT sample stored temporarily in, for 
example, RAM. 
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It will be noted that the frequency domain sequence spectrum S k 262 is 
essentially flat with slight perturbations about a nominally constant magnitude. 



The sequence spectrum S k 262 is compared with the channel sequence 
spectrum Y k 254 in a cross-correlation function 266, such as realised by Wiener 
filter frequency response, namely 

K 

2 2 , to produce an impulse response spectrum H k 268 that can 
+ <7 



be stored in memory for the duration of at least the related slot and longer if 
desired. Clearly, by storing the impulse response spectrum H k 268 for an 
extended period, processing overhead is reduced at the expense of accuracy in 
any correction of the variations in the channel due to Doppler. As will be 
understood, time alignment of the received and transmitted training sequences 
is required prior to cross-correlation. 



Optionally, the impulse response spectrum H k 268 may be further refined and so 
it is subjected to an IFFT function 270 to generate an estimated time domain 
representation h k thereof having an arbitrary number of taps/channels that can 
be used as a matched filter, if desired. In other words, each tap is representative 
of a weight for the channel impulse response. Within the time domain, trimming 
of the number of taps reduces the channel impulse response to manageable 
processing levels, with the reduction reflecting the greatest expected channel 
impulse response (CIR) duration (presently understood to be seventy-one chips 
in length in 3GPP, although this number is an arbitrarily selected cut-off). From 
the reduced number of taps in the channel impulse response, most values will 
approximate to zero and so a control processor 272, based on a level or number 
threshold, selects only those N significant taps that most noticeably reflect the 
channel to produce a best estimate channel impulse response h' k 274. Selection 
of the N most significant taps may simply force lesser taps to be forced to zero. 
By limiting the channel impulse response to the most significant taps, processing 
is reduced and a final error in the channel limited. In the latter respect, it will be 



NORTEL CONFIDENTIAL PROPRIETARY 



FPUS00670 



-21- 



understood that all tap values will have an associated error from thermal noise, 
channel interference and other adverse channel and processing artefacts, so the 
compounding of errors is mitigated. 

5 A complete understanding of the principles of the present invention will now be 
gleaned from the following description of specific signal processing interactions 
and transmission techniques. 

Concerning transmission of a training sequence according to a preferred 

10 embodiment of the present invention, reference is made to FIG. 4. For the sake 
of explanation, the structuring of base station transmissions (in the timing frame 
of FIG. 4) is assumed to reflect a three site system in which each of three base 
stations includes a single transmit element. Of course, the basic principle 
employed in the specific embodiment of FIG. 4 can be developed for multi- 

15 element transmit arrays from two or more base stations. Essentially, FIG. 4 is 
representative of successive training sequence transmission bursts from each of 
three base stations (Base TX1, Base TX2 and Base TX3). In generality, the 
number of training sequences transmitted from all transmit elements follow a 
squared law, i.e. for three base stations each having single radiating elements 

20 then there would be three successive bursts (ti-t 3 ) separated, preferably, by 
guard periods (gi and g2) or cyclic prefixes. Bursts from the various base 
stations are preferably time aligned at transmission, although a windowing 
function at the receiver (that steps through a chip-wise time correlation) may be 
employed to time align the transmitted sequences s n k. Since channel delay is 

25 generally no greater than about 20p.s, then the receiver need only step through a 
chip count commensurate with the maximum channel delay. For example, at an 
operating frequency (chipping rate) of one megahertz (1MHz), 20|is correspond 
to twenty chips. 

30 During each burst (ti-ts), each base station transmits a finite (predetermined) 
length of training sequence s, which training sequence may vary from 1-chip or 
more. Typically, each burst will be tens of chips in length and preferably at least 
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about sixty chips in length and most preferably in excess of about one hundred 
chips. Clearly, the length of the training sequence is restricted to manageable 
levels to avoid inefficient use of communication resource for more important data 
transmissions. The nomenclature adopted in FIG. 4 with respect to the bursts 

5 follows the convention that the first subscript (n) represents an identity of the 
transmitting base station whereas the second subscript (k) represents the 
sequence number. In the context of Base TX1, FIG. 4 is illustrative of three 
sequences sn, S12 and s 13 during bursts ti, t 2 and t 3 , respectively. A channel 
sequence spectrum Y k is formed (for each burst (ti-t 3 )) from the summation of 

10 the individual sequences sn, S21, s 31 from the base stations transmitting in that 
burst. 



For channel estimation, it will be appreciated that the number of samples 
required is dependent upon bandwidth, e.g. with a symbol rate of 284kbps (kilo 
1 5 bits per second), GSM systems require about five samples of the channel. 



It is well known that DFT's have a wrap-round effect that is equivalent to the data 
being cyclic. Consequently, in (contiguous) observation windows each encoded 
with the same PN sequence, DFT's can be used to perform a discrete 

20 convolution operation. Circulant matrices are equivalent to these DFT 
operations, which circulant matrices can be defined as follows. In a system in 
which the training sequence is S = {s 0j $2, $2, s 3 , s *h the cyclic extension is the 
channel impulse response is H = {ft 0 , /?i , h 2 , /? 3l h 4 } and the observation window 
arranged to hold data Y = {yo, yi, y2, y3, y^, then Yis the cyclic convolution of S 

25 with H: 

N 

77=0 

with N=4. If the associated DFT's are the sequence spectrum Sk, the channel 
impulse response H and the channel sequence spectrum Yk is then: 

Y k =S k ®H k 

30 where ® is a point-by point-product. 
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In developing a solution to the problem of channel estimation in the presence of 
interference or multi-user channel estimation (in either non-dispersive or highly 
dispersive channels), the inventor has appreciated that the fundamental 
5 problems associated with channel isolation/resolution can be expressed 
algebraically in similar terms to multi-user detection algorithms, as referenced 
above by J. Blantz et al and A. Klein et al. 

In the Steiner code method described above, to estimate M channels the length 
10 of the PN code must be at least M times the maximum channel dispersion. In an 
N base station situation, each base station with M elements, while the Steiner 
method can continue to be used, the PN code length must be at least MN times 
the channel impulse response duration and each element of each base station 
would transmit an end-around rotated version of the code with a different unique 
15 offset, in the multiple base station situation, the differing round-trip propagation 
times to the subscriber unit must also be taken into account in the channel 
model so codes becomes increasingly extended. A more flexible option for 
allowing expansion of the number of serving base stations is to use an 
alternating sequence of different Steiner codes, each of which on its own is only 
20 just long enough for the number of elements in a single base station to be 
resolved. 

To appreciate the present invention, it is preferable to consider a communication 
system having N active/serving base stations for a call, with each base station 

25 having a single transmitting element. If the N base stations each with a single 
element are received by a single element antenna of a subscriber unit, such as 
cellphone or terminal, then if the N bases simultaneously send N strings of N 
different pseudorandom bursts and such that there are N 2 different sequences 
altogether, then the individual channel impulse responses can be resolved at the 

30 terminal. In practice, the present invention is particularly interested in a three- 
base station situation which, statistically, occurs most frequently. 
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The base stations are preferably synchronised to a time accuracy much less 
than the duration of any cyclic prefix (which is typically of the order of about 10^xs 
to 20ps for urban microwave propagation). The m Xh received burst is now the 
convolution: 

where b is the thermal noise component and K is the maximum dispersion in 
chips. After transformation into the DFT (frequency) domain, the equivalent, 
numerically efficient function can be expressed in the point-by-point product 
form: 

10 7 (w) = f^S { „ m) ®H n +B {ni) 

where 

, m \DFT 

15 

/ \ DfT 

Thus, rather than performing a convolution of the channel impulse response h n 
with the m th sequence s mn from the n ih base transmitter element, the present 
invention applies a DFT operation to both the received sequence y n and the 

20 various local replicas of the training sequences s mm and then does a 
multiplication together of the corresponding frequency bins of the received signal 
sequence Y n (being the signal sequence s n as received through a channel) with 
the corresponding frequency bins of the replica of the training sequence S mn . 
Summation (i.e. concatenation) of the individual multiplication products 

25 associated with each channel provide a composite frequency response for the 
channel, the composite frequency response allowing, in the time domain, 
generation of the channel impulse response for the channel. The process of the 
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present invention generally involves periodic buffering, for example, of frequency 
bin data and vector results pertaining to channel impulse responses, as will be 
appreciated by the skilled addressee. 

5 In other words, taking an exemplary 128-point Fourier transform of both the 
received sequence y with the training sequence s yields one hundred and 
twenty-eight corresponding frequency bins in the frequency domain. 
Multiplication together of corresponding frequency bins of the received sequence 
with the same frequency bin of the signal (training) sequence s generates a 
10 matrix containing solvable linear equations in H and S. Summation (i.e. 
concatenation) of the individual multiplication products associated with each 
channel then provides a composite frequency response for that channel, which 
composite frequency response can be used in the time domain (through an 
inverse FFT) to provide a channel impulse response. 

15 

It will be appreciated that the channel is assumed constant for the duration of the 
burst events with respect to slow moving subscriber units, e.g. mobile units 
having a velocity of less than about ten kilometres per hour (10kph). 

20 In the DFT domain, the /c* frequency bin received for the first burst can be 
expressed as the sum of components from the base stations and their channels: 

Similarly, the k* h frequency bin for the second burst can be expressed as: 

7 2 w =SW > 

25 And, finally, the fc" 1 frequency bin for the third burst is: 

Rearranging the tf h frequency bins provides a set of solvable linear equations in 
terms of H and S, namely: 
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which can be written in the matrix vector form at each frequency k bin (which for 
clarity is omitted): 

5 It is now possible to define matrices and vectors containing all the transmitter 

elements and bursts: 
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These equations are solved using a minimum mean square error (MMSE), e.g. 
1 0 Weiner filter, method providing: 

H ^{s H S + Cov{B)) X S H Y 

Here Cov(B) is the covariance matrix of the thermal noise, which is often 
approximated by an identity matrix a z l. The equation immediately above can be 
considered as two entities, namely as a matched filter S H Y arranged to minimize 
15 the signal to noise ratio but not to remove mutual interference between channel 
estimates, followed a decoupling matric {S H S+Cov(B))~ 1 that removes the 
coupling between different channels caused by the non-ideal nature of the PN 
sequences. The inverse DFTs for each radio channel are therefore of the form: 





iDFT 
=> 
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1, 00 
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20 where z is the frequency bin (or its time domain equivalent). 
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The decoupling matrix term S H S is a matrix product, where S H is the conjugate 
transpose of the sequence, i.e. the transposed matrix with all elements complex 
conjugated. The covariance matrix <SS H > is the time averaged outer product of 
the noise vector, or in element-by-element terms: 



<BB h 



r <b x b[> <b x b* 2 > <b { b* 3 > 
<b 2 b* > <b 2 b* 2 > 
<b 3 b* > <b?b* 2 > 



In terms of each channel concatenation of the numerous frequency bins allows 
generation of a time domain vectors of the channel impulse response for a 
channel, the time domain vectors generated from an inverse DFT process within 
10 a suitable processing block, such as an ASIC or control processor 150 (of FIG. 
1). 



It is now possible to extend the basic concept of multiple channel estimation to a 
situation in which each base station has a plurality of transmit element; this 

15 scenario arises in, for example, systems employing space time coding (STC). It 
will be appreciated that STC is a collection of techniques for transmitting a 
number of parallel data streams from multiple transmit antennas with the aim of 
increasing bit rate for a given total transmitter power. In STC, it will be 
appreciated that the 2x2 frame structure is expanded to an N-element array, 

20 usually having at least four transmit elements. For example, in a four element 
STC transmitting antenna we might organise the user data in blocks of four 
symbols and transmit them in different order from the four elements via 16 
channels to 4 receiving elements. This is a natural extension of the two element 
STTD procedure, the main difference from STTD being that multiple receiving 

25 elements are essential to get a bandwidth expansion in STC whereas STTD is 
aimed purely at improving the reliability of the received data and is not 
concerned with bit rate expansion. Thus STTD will function with a single receive 
antenna, while STC requires a minimum of two. The principle difference between 
an STTD arrangement (that is discussed in detail below) and STC is that each 
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component of the block vectors and block matrices (for STTD given below) 
contains an increasing number of elements (i.e. CIR samples, transmitted data 
samples and received data samples) consistent with the adopted coding method 
in each STC derivative. 

5 

Referring briefly to FIG. 5, which is generally complementary to FIGs. 1, 3 and 4, 
there is shown a schematic block diagram of the channel estimation process and 
architecture 500 of a preferred embodiment of the present invention, the process 
based on a two base station (BS1 and BS2) transmission environment having 

10 two training bursts. The training sequences s x (co) are, in this exemplary case, 
Gold codes, with individual first 502 and second 504 training sequences 
transmitted in a first burst 506. The first training sequence 502 is subjected to 
first channel 508, whereas the second training sequence is subjected to a 
second channel 510. Both channels are combined at the receiver 512 and are 

15 further subject to a first noise component (N^(<o)) 513. In this exemplary case, 
first and second matched filters 514 and 516 provide a time reversed conjugated 
inversion of the respective first and second training sequences 502 and 504, 
with outputs from the first and second matched filters 514 and 516 provided to a 
frequency domain coupler 518 supporting the frequency domain decoupling 

20 matrix (S H S+Conv(B))' 1 . Third 520 and fourth 522 training sequences transmitted 
in a second burst 524 are respectively subjected to first channel 508 and second 
channel 510. Both channels are combined at the receiver 512 and are further 
subject to a second noise component 526 (typically different to the first noise 
component). In this exemplary case, third 528 and fourth 530 matched filters 

25 provide a time reversed conjugated inversion of the respective third and fourth 
training sequences 520 and 522, with outputs from the third and fourth matched 
filters 528 and 530 provided to the frequency domain coupler 518. The 
frequency domain coupler then provides channel impulse response vectors 
Hi(co) and #i(g>) (reference numerals 532 and 534) for the first and second 

30 channels, respectively, based on the foregoing point-by-point manipulation of 
individual (frequency domain) frequency bins. 
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The preferred embodiment of the present invention therefore shows how 
accurate channel estimation is possible through the use of Wiener frequency 
domain MMSE deconvolution combined with frequency domain spatial 
decoupling matrices, with quasi-orthogonal Gold codes allocated to base 
5 stations and their antenna elements. 

The principles of the present invention can also be applied to resolve multiple 
channels by a single training burst; this has particular application in OFDM 
systems where it may not be possible to transmit Steiner codes because of 

10 hardware constraints. The present invention can be modified, in the fashion now 
shown, to operate with a single arbitrary training sequence yet still resolve 
multiple channels in a least squares sense. A single training burst situation is 
algebraically undetermined in the two channel case. If the burst is of length L } 
there are 2L unknowns in the two channels which cannot be determined from L 

15 observations at one antenna. However, if there are also VkL linear constraints 
associated with each channel impulse response the equations become solvable 
again. The simplest constraints are that the channel impulse responses {h 0 ...hk} 
are of finite duration. If K< YzL then we write the equations as follows. Let S m be 
the circulant matrices representing the training bursts and H m be the IR vectors 

20 for channels #1 and #2. However, with a matrix solution not be circulant, the 
ability to use FFT to solve the following equation is lost, but resolution of the 
channels is still possible with use of a matrix of the form: 

S \\ ^10 ^21 ^20 
Y— S n S u S 22 #21 

0 s n 0 s 22 
0 0 0 0 

In generality, with respect to two parallel channels using a single training burst, if 
25 there are a large number of linear constraints C attached to the channel impulse 
responses (such as finite durations) then a number of channels can be solved 
with a single training burst. Algebraically, if the there are N c constraints on each 
channel vector of L elements and C channels, the sum of the number of 
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independent degrees of freedom is C(L - Afc). If the length of the training 
sequence is N T > C(L - Nc), then all of the independent degrees of freedom can 
be solved; this method is used in Steiner codes to estimate multiple channels 
using a single training sequence 

5 

Whilst benefits are attained through use of the channel impulse response 
estimation scheme described so far (in which the number of processing blocks is 
equal to the total number of channels to avoid code contention problems), some 
processing efficiency (i.e. DSP load) associated with Steiner codes is lost since 
10 there is a requirement to invert a matrix in each frequency bin whose dimension 
is equal to the number of channels, i.e. the number of transmitting cells times the 
number of base antenna elements per cell. 

In a further preferred embodiment of the present invention, therefore, the training 

15 sequence transmissions of FIG. 4 are implemented as Steiner cyclic training 
sequence variant; this is shown in FIG. 6 in which successive bursts 602 and 
604 from different base stations BS1 and BS2 contain repeated sequence 
portions of an entire training sequence. This particular embodiment of the 
present invention offers an optimum combination of frequency domain Wiener 

20 filtering and MMSE spatial decoupling for the estimation process. The high 
accuracy of channel estimation that the combination of techniques offers is a 
desirable if not essential for coherent null steering and soft hand off in CDMA 
and OFDM. The sequences {sn}, {S21}, {Si2>, {S22} are four linearly independent 
random sequences transmitted in pairs from the two antennas as shown. These 

25 sequences transform into the DFT domain as Sn, S21, Si2> S11. Similarly the 
channel impulse responses are Hi, H 2 in the DFT domain (for the two channels 
in the exemplary case of a two base station system having single transmit 
elements). The receiver (e.g. receiver chain samples the two received blocks 
and, via a DFT function, transforms to the Frequency domain as Y1 and Y 2 . The 

30 required guard intervals, which prevent dispersion contamination by unknown 
data, are already provided in the Steiner sequences and additional, i.e. 
separate, provisioning of guard periods is not therefore required. The two 
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received blocks are superpositions of convolutions of the transmitted sequences 
with the channels and are modeled by point by point equations in the DFT 
domain as follows: 
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5 Since the S matrix has no particular orthogonality properties, matched filtering is 
not advisable. Instead, if the following MMSE matrix operation is performed at 
each frequency; the result combines spatial decoupling with a Wiener frequency 
weighting function and overall gives an MMSE solution. It requires a 2x2 matrix 
inversion in each frequency bin for two cell operation (with each base station 
1 0 having a single transmit element): 
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By combining Steiner cyclic pilot rotation for transmissions from the transmit 
elements of each base station according to the present invention, 
interference-free mobiles near a centre of the cell are allowed to resolve their 
15 channels in the usual (conventional) way. At the cell periphery where 
interference becomes significant, the training sequences need to be expanded 
to provide enough degrees of freedom such that the simultaneous equations for 
all of the channels become algebraically determinate; a repeating pattern of the 
Steiner PN sequences allows this to occur. 

20 

In the exemplary case where a receiver observes six sequences, as shown in 
FIG. 7, the result in the Fourier domain is that there are sufficient degrees of 
freedom to solve for six channels to {h G }. It should be note that, in FIG. 7, the 
introduction of sign reversals allows for distinction to be drawn between codes in 
25 successive bursts, with the mechanism of sign reversal therefore reducing the 
overall number of Steiner codes required to implement a system according to the 
present invention. The equations relating the received signal Y to the sequences 
U, V, W and the channels H in the frequency domain are: 
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However this is an inefficient algebraically redundant representation that can be 
simplified as follows. At the mobile terminal, for example, the signal from a given 
antenna at a base station is the circular convolution of the PN sequence {s} with 
5 the channel impulse response {h}: 

N-l 

y k = Z"V;mod A- * = 0- . .2" -1 

The cyclic shift of the PN sequence for second base element in the Steiner 
method can be replaced by a reverse shift of the channel impulse response to 
give exactly the same signal at the subscriber unit: 

N-l 

Vic ~ y ^d S k-j-%MODN^ l j 

10 

AM 

= S S k-jhj+% MOD N 
7=0 

Cyclic shifting of channel impulse responses /? 2 , h A and he (of FIG. 7) takes into 
account the reverse shift as detailed immediately above with the effect that the 
six sequences received at the antenna of the subscriber (or other addressed 
unit) are Y 2 , Yb, Y 4i -Y$ and -Y 6 . Moreover, the six sequences are obtained 

15 from a repetition of only three Steiner codes that are individually and uniquely 
assigned to each base station, i.e. both transmit elements send identical 
sequences in the same burst, with the subscriber unit generating a cyclically 
shifted channel impulse response ft' x for one of the two radiating elements of the 
base station's transmit array. This reverse shift manipulation results in only three 

20 different PN sequences being utilised in this particular embodiment, with the 
channel impulse responses of FIG. 7 taking the form {fr\ + // 2 ), Qte + tit) and {h$ 
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10 



15 



+ /?* 6 ) to provide the received sequence Y\, Y 2 , Yb> Y 4 , -Y 5 and -Ye- The 
determinant of the sign change matrix of this solution is 4, as shown below: 
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with Fourier transforms: 



yi=> Y i y 2 =^Y : 

hj^Hj h 2 =>H 2 
u=>U v=> V w 
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then a reduced set of more easily solvable equations in the DFT domain is 
obtained, namely: 
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Provided a resultant matrix in this set of equations is well conditioned and 
invertible, then the set of channel impulse responses can be solved, A sampling 
design for the repetitions that offers near-orthogonality between repeated 
sequences is, however, required, such that different bases can be resolved 
without ill-conditioning of the equations. Such designs, as used for space-time 
codes, are readily available in statistics literature. 



By employing the preferred embodiment of the present invention where Steiner 
codes supplement to point-by-point DFT multiplication mechanism, accuracy in 
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estimated channel impulse responses is limited by channel noise rather than 
imperfection in the codes. 

The preferred embodiment that uses Steiner codes for m base antenna 
5 elements can be used in multiple cells and MMSE channel estimation is 
obtainable for a number n of bases giving a total of nm downlink channels. 
Different cells will use a different basic PN code, whereas the elements within a 
single cell will use the Steiner cyclic offsets of a single PN code. The extra 
workload (processing overhead) involved is the inversion of an nxn matrix in 
10 each frequency bin. The inclusion of the use of Steiner codes to supplement the 
basic inventive concept of Wiener frequency domain MMSE deconvolution 
combined with frequency domain spatial decoupling matrices results in the 
possibility of allocating only a single training sequence to each base station 
provided that the training sequence is of sufficient length to encompass all the 
15 multiple time-translated channel impulse responses. 

Embodiments of the invention may be implemented as computer program code 
encoded on a computer program product for use with a computer system. It is 
expected that such a computer program product may be distributed as a 

20 removable medium with accompanying printed or electronic documentation (e.g. 
shrink-wrapped software), preloaded with a computer system or distributed from 
a server or electronic bulletin board over a network (e.g. the Internet or World 
Wide Web). A series of computer instructions can therefore either be fixed on a 
tangible medium or fixed in a computer data signal embodied in a carrier wave 

25 that is transmittable to a computer system using wireline or wireless 
transmission techniques. For example, the present invention can find application 
as downloadable code in a software radio. The removable (i.e. tangible) medium 
may be a computer readable media, such as a diskette, CD-ROM, DVD-ROM or 
RAM, fixed disk, magneto-optical disks, ROMs, flash memory or magnetic or 

30 optical cards. The series of computer instructions embodies all or part of the 
functionality previously described herein with respect to the system. 



NORTEL CONFIDENTIAL PROPRIETARY 



FPUS00670 

-35- 

Software embodiments of the invention may be implemented in any conventional 
computer programming language. For example, preferred embodiments may be 
implemented in a procedural programming language (e.g. "C") or an object 
oriented programming language (e.g. "C++"). 

5 

Although the preferred operating method is realised by general or specific- 
purpose processor or logic circuits programmed with suitable machine- 
executable instructions, hardware components may possibly be used to 
implement certain features of the present invention. Of course, the present 
10 invention may be performed by a combination of hardware and software. 

It will, of course, be appreciated that the above description has been given by 
way of example only and that modifications in detail may be made within the 
scope of the present invention. For example, while the preferred embodiment is 

15 has been described in relation to a multi-cast system in which each base station 
has a single transmit element, the present invention is clearly applicable to multi- 
element transmitters and to single isolated cells supporting transmissions from a 
single base station having a transmit array. Furthermore, while the preferred 
embodiment has been described in the context of a CDMA system, the present 

20 invention is generally applicable to any dispersive or imperfect channel that acts 
to corrupt symbols, including systems employing TDD and FDM techniques and 
other types of links, such as digital subscriber lines (DSL) and all digital 
transmission systems. 

25 Indeed, although the technique of the preferred embodiment of the present 
invention has been described for a CDMA application where the number of chips 
in the training burst is easily identified, the principles of the present invention 
apply to OFDM, wavelets (e.g. Daubechies's wavelets) and any similar 
orthogonal or near-orthogonal signaling waveforms. With specific regard to an 

30 OFDM signal, such an OFDM signal is treated as a single wide-band modulated 
carrier in the time domain and an estimate of the time domain channel impulse 
response is made in an identical way to that described in the preferred 
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embodiment. If a constraint is made that the duration of the impulse responses is 
small, the equations become algebraically soluble (or "observable" in the context 
of control theory jargon) and in this way a large number of parallel channels can 
be estimated using a single training burst. 

If the duration of the OFDM training burst is N s sub-carrier symbols and there 
are N F sub-carriers then the burst effectively has A/ S *A/ F samples or chips which 
is quite a substantial number. In practice, this means that multiple channel 
estimation in OFDM might not be too difficult. The lower limit for the duration of 
the training sequences would be that it must be capable of supporting an 
orthogonal, or at least a linearly independent, set of codes for say m base 
station antennas seen from the terminal (in one or more cells). Superficially, this 
would require m full sub-carrier slots and for HSD with 1 50 (is slots we need a 
duration of 150m \is. However the wideband methods described here only 
require that the duration of the training sequence is > m times the duration of the 
channel impulse response so, for a channel with t jas dispersion and signal 
bandwidth B MHz, the duration of the training sequence need only be xBm 
samples. For example, if the signal bandwidth is 6 MHz, the channel dispersion 
is 5|its and there are seven base station antennas with channels to be 
determined, the minimum duration of the training sequence is 210 samples or 
6x5x7/6 = 35^is, i.e. only 1/5 of a nominal sub-carrier symbol. In practice, 
however, this duration would need to be increased significantly (and preferably 
about ten-fold or so) to improve estimation accuracy. 

In transmission bursts from base stations having single transmit elements, cyclic 
prefixes could be replaced by blank (zero) carriers to provide isolation between 
burst sequences. 

It should be understood that reference to a base station refers to serving 
infrastructure irrespective of the form of air-interface adopted by the exemplary 
CDMA-based architecture of the preferred embodiment. In other words, the term 



NORTEL CONFIDENTIAL PROPRIETARY 



FPUS00670 



-37- 



base station should be construed broadly and, where necessary, to include 
references to communication relay devices. While the main applications are in 
channel estimation on a downlink for CDMA and OFDM in order to apply 
coherent soft hand-off and interference cancellation in locations where two or 
more base stations are at similar ranges and the weaker base stations, if 
ignored, cause mutual interference, the preferred embodiments are not intended 
to be limiting to these specific arrangements. 



We claim: 
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